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The theory of specular X-ray reflectivity from a rough interface based upon the reflection function 
method (RFM) is proposed. The RFM transforms the second order differential equation for the 
wave amplitude into the non-linear first order differential equation of Riccati type for the reflection 
function. This equation is solved in the approximation of the abruptly changing potential, which is 
justified for the typical angles of X-ray reflectometry. The reflectivity is represented as a series. The 
first term of this series reproduces the Nevot-Croce approximation and second one gives the phase 
correction for greater angles. It is shown that the phase correction can be used to obtain the degree 
of interface asymmetry. The X-ray reflectometry model profiles for Fe/Cr superlattice are used to 
illustrate the method. 



The X-ray reflectometry is a useful tool for studying 
surface and interface structure in thin films and multi- 
layers. Usually, the rough surface X-ray reflection is an- 
alyzed in the framework of the plane- wave Born approxi- 
mation (PWBA) or the distorted-wave Born approxima- 
tion (DWB A) Jl|, H . In this work, we apply the reflection 
function method (RFM) j| to the specular X-ray reflec- 
tion from a rough surface or interface. 

Let us consider the X-ray reflection on a non-ideal in- 
terface structure. We assume that this structure is homo- 
geneous along the surface which is parallel to the (x, y) 
plane and the media can be characterized by its dielectric 
susceptibility x( z ) depending only on the normal coordi- 
nate z, where x( z ) ~ y when z — ► ±oo. The change of 
the material occurs only in the z-direction perpendicular 
to the surface. Then one has to solve the one - dimen- 
sional Helmholtz equation 



electric field E(z) in the form 



+ k 2 sin 2 9 ) E(z) + k 2 x(z)E(z) = (1) 



dz 2 



Here E(z) is the electric field in the medium, 6 is the 
incident angle and k — 2tt/X, A being a wave length of 
radiation. As the first step, we need to evaluate the scat- 
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interface between two subsequent layers, which are de- 
noted as 1 and 2. The RFM starts from the transforma- 
tion of the linear second order differential equation (|l|) 
for the wave amplitude E(z) into a non- linear first or- 
der equation of Riccati type for the reflection function 
B(z). This transformation is not unique and can be per- 
formed in a number of different ways. An advantage of 
the RFM is that the perturbation expansion carried out 
in the framework of this scheme gives more rapid conver- 
gence in comparison with the conventional Born series. 
In particular, the first order approximation easily enables 
one to go beyond DWB A. 



We denote q(z) — 2k \ sin 6 + x( z ), an< 4 represent the 



E(z)=q- 1 / 2 (z) 



A{z) exp ( — J q(x)dx ) + 



C(z) exp 



q(x)dx 



(2) 



where A[z) and C[z) are amplitude functions. In addi- 
tion, we apply the following condition: 



dz 



E(z) 



= V/2 



A(z) exp 



-C(z) exp 



q(x)dx 



q(x)dx 



(3) 



The reflection function B(z) is defined as B(z) = 
C{z)/A{z). Taking into account the continuity of E{z) 
and Eq. (0) one can prove that B(z) satisfies the first 
order nonlinear differential equation 



dz ^ ^ %q(z) 



exp [ 1 J q{x)dx 



£? 2 (z)exp^— i J q{x)da 



(4) 



Eqs. (|l|) and (m) should be supplemented by the bound- 
ary conditions. For example, the choice of B(+oo) = 
corresponds to the X-ray beam, incident from z < 0, and 
in this case a reflection coefficient r\\ is given by rela- 
tion m = B(— 00). We also introduce into consideration 
the dimcnsionlcss functions g±(z), which are related to 
XO) via equality x( z ) ~ X± = ±(X- - X+)g±( z )- The 
function g~(z) — > 0, when z — > —00, and g-{z) — > 1, 
if z — > +00 (See Fig.l). The functions g±(z) obey the 
relation g+(z) + g~{z) = 1. One can regard g±{z) as 
a "shape" of the interface, which reproduces the grad- 
ual transition from the first layer to the second one. We 
shall call interface "symmetric", if (-^-)g-(z) is an even 
function of z, otherwise interface is "asymmetric". 

In case of grazing incidence angles Eq. (^) can be solved 
in the approximation of the abruptly changing potential. 
The small parameter e of this expansion is defined as e = 



2 



aq c /2ir, where a is characteristic length corresponding to 
the variation of the potential and q c = max\q(z)\. In 
the X-ray reflectometry studies a is of the order of mean- 
root-square interfacial roughness a = 2 — 8 A(See 0) and 
q c ~ (47r/A)sin0. Therefore the condition e < 1 holds 
over the scattering angle region (0 < 8 < 4°). These 
estimates make it possible to find the solution of Eq. (3) 
in the form 



where 



B{z)=B (z)exp(p(z)), 



B {z) = (q(z)-q 2 )/(q(z)+q 2 ), 



(5) 



92(i) = 2fc y sin 9 + x±, 



and 



/?(z)=£/? n (z)e" 

ra=l 

The function Bq(z) corresponds to the boundary con- 
dition B(+oo) = and it gives the Fresnel reflection 
coefficient fn = (qi — q2)/{qi + 92) from an ideal sharp 
interface. The series (3{z) yields the corrections due to 
the interfacial non-ideality. 

The use of ansatz (||) is the essential step in the deriva- 
tion. It enables us partially to sum up the reducible parts 
of the expansion B{z) in powers of e, so that the coef- 
ficients P n (z) are associated with the irreducible terms 
only. The series (3(z) can be found by means of subse- 
quent iterations from Eq. (Q). It turns out that, at each 
step n, one encounters the only linear inhomogeneous dif- 
ferential equation for (3 n (z). The details of this derivation 
will be presented elsewhere. As a result, up to the third 
order of e, the elements of the matrix S12 can be written 
in the form 



r| 5 <<xi> ( iqxS - iqi^c 2 



T22 



(0) 



r 22 ex P ( ~^2<5 - -?1 <?2C 

n2 [-(qi + 3g?)^! - (ii - ql)3\ ^) (?) 



*12(21) = ^12(21) CX P^2*(' ?1_92 ) (5+ 8^ 1 " 92 ) ° ( 8 ) 

+ \i{qi - qi) 2 [(91 - ©Om? + (qi + 92)^2] v 3 ) 

Here r F , t F are Fresnel's reflection and transmission am- 
plitudes and parameters 5, <r, and jJL\{2) are expressed via 
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FIG. 1: The linear segment form of the profile g~(z), corre- 
sponding to the interface of a width 2a. The "symmetric" 
case is shown by the dashed line, and the "asymmetric" one 
is depicted by the solid line. 



g±(z) as follows 
S = 



z—g-(z)dz 
az 



(9) 



Ml(2) 



/+OC /- + 0O 

9-{z 1 )dz 1 I g + (z 2 )dz 2 (10) 
-oo J Z\ 

/ (3) (-,+,+)T/ (3) (-,-,+)l/4^ 3 (11) 



where 

I (3) (-,±,+) = 



+ 00 Z' + OC 

g-(zx)dzi / g±(z 2 )dz 2 

— oo J Z\ 

g + (z 3 )dz 3 
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Consider now the physical meaning of Eqs. (^j|). First 
of all, the phase shift 5 arises due to transmitting electro- 
magnetic wave in the non-uniform interface region. This 
phase shift is equivalent to "effective" increasing in the 
thickness of layer 1 to value 8: z = z' — S. (See Fig.l) 
In the process of the numerical treatment of the X-ray 
reflectometry profiles this fact enables one to adjust the 
ratio between the layers' thickness' in the periodical cell 
of the superlattice in order to obtain the best fit to ex- 
perimental data. The second order correction to the am- 
plitudes r F , t in Eq. (6-8) reproduces the well-known 
Nevot-Crocc [^| approximation. The magnitude er has 
the meaning of the root-mean-square interfacial rough- 
ness and it is given by Eq.(|To|). 

The phase correction corresponding to the third or- 
der terms of the expansion is a new feature in the ques- 
tion. In addition to a it contains two extra parameters 
Hi and \X2- We found, that \i 2 is in general nonzero for 
a wide set of profiles x( z ) whereas \x\ does not vanish 
in case of the asymmetric interfaces only. Hence, this 
property may be used to define \i\ as the measure of the 
interface asymmetry. The parameter [12 has no such an 
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FIG. 2: Model X-ray reflectivity profiles for multilayer struc- 
ture Al 2 O 3 /Cr(70A)/[Fe(20A)/Cr(9A)] 8 calculated without 
asymmetric phase corrections (points), and with asymmet- 
ric phase correction (fif — 0.2, solid line). Wave length 
A = 1.789A. 

evident meaning as fix- But we may note, that if nec- 
essary, it can be eliminated from the consideration. It 
follows directly from the structure of Exp. (6-8). Since 
the factor — g| = 16(7r 2 /A 2 )(x- — X+) does not de- 
pend on 6, the terms, proportional to /xfj give mere 
the additive contribution to S. One may take it into 
account by omiting the terms with and substituting 
<5^ + 16( 7 r 2 /A 2 )(x--X + Vi 

To test the obtained approximation we exploited 



the symmetric Epstein profile g^(z) = (1 + e z / a ) 1 
for which the exact solution is known. In this case we 
obtained a = {n/VS)a and [i\ = (3V3/2tt 3 )C(3) « 0.100. 
Assuming further /ijj = 0.1 the model X-ray profile cor- 
responding to the Al 2 O 3 /Cr(70A)/[Fe(20A)/Cr(9A)] 8 
multilayer have been calculated, taking into account the 
possible asymmetry /ii in the interfacial structure. Pro- 
vided the matrices Sk,k+i are known, the solution of the 
Eq. (|l|) and, hence, the scattering matrix S of the whole 
multilayer is found by means of recurrent scheme [p|. 
The results obtained are shown in Fig 2. In agreement 
with Eqs. 1|) the phase correction becomes essential 
with the increase of the incident angle and it provides a 
more adequate description of the reflectometry spectrum 
for the scattering vectors in the range from the first to 
the second Braggs' peaks. The another sample profile is 
shown in Figs. 3 and 4. It corresponds to the structure 
Cr203(<T=3A)/Fe(a=2A)/[Cr(a=2A)/Fe(o-=2A)]n/Cr 
(<7=3A)/Al 2 3 ((7=lA); Cr 2 O 3 (m=0.77)/Fe(//i=1.25)/ 
[Cr(// 1 =-0.0045)/Fe(/ii=-1.0)]n/Cr(/i 1 =-0.56)/Al 2 O 3 
(/ii=2.4). The more exhaustive account and the details 
of our numerical algorithm will be presented elsewhere. 
We would like to emphasize that the form of the scatter- 
ing matrix as given in Eqs. is rather general, i.e., it 
is irrelevant to the precise form of a reflectivity profile. 
Thus it provides the unification description of a large 
variety of possible symmetric as well as asymmetric 
interfaces. 

Summing up, we have developed the theory of specu- 
lar X-ray reflectivity from a rough interface based upon 
the reflection function method. By using the approxima- 
tion of the abruptly changing potential we have found 
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FIG. 3: Model X-ray reflectivity profiles for multilayer struc- 
ture Al 2 O 3 /Cr(70A)/[Fe(70A)/Cr(9A)]i2 calculated without 
asymmetric phase corrections (sig), and with asymmetric 
phase correction (phase). Wave length A = 2.070A. 
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FIG. 4: Model X-ray reflectivity profiles for multilayer struc- 
ture Al 2 O 3 /Cr(70A)/[Fe(70A)/Cr(9A)]i2 calculated without 
asymmetric phase corrections (sig), and with asymmetric 
phase correction (phase). Wave length A = 2.070A. 
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the phase correction to the reflectivity due to interface 
roughness and asymmetry, which is essential for the de- 
scription of the X-ray reflectivity spectra for greater in- 
cident angles. 
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